CLNS 93/1244 
LBL-34573 
hep-ph/9308364 



Radiatively-Induced First-Order Phase Transitions: 
The Necessity of the Renormahzation Group 

Mark Alford 
Laboratory of Nuclear Studies, 
Cornell University, 
Ithaca, NY 14853 

and 

John March-Russell 
Theoretical Physics Group, 
Lawrence Berkeley Laboratory, 
1 Cyclotron Rd, 
Berkeley, CA 94720 

Abstract 

We advocate a (Wilson) renormalization-group (RG) treatment of finite- 
temperature first-order phase transitions, in particular those driven by 
radiative corrections such as occur in the standard model, and other 
spontaneously-broken gauge theories. We introduce the scale- dependent 
coarse-grained free energy S'a[0] which we explicitly calculate, using the 
Wilson RG and a (4 — £)-expansion, for a scalar toy model that shares 
many features of the gauged case. As argued by Langer and others, the 
dynamics of the phase transition are described by 5'a[0] with A of order 
the bubble wall thickness, and not by the usual (RG-improved) finite- 
temperature effective action which is reproduced by S'a[0] for A — > 0. We 
argue that for weakly first-order transitions (such as that in the standard 
model) the (4— £:)-expansion is necessary to control an inevitable growth of the 
effective scale-dependent coupling towards the strong-coupling regime, and 
that diagrammatic resummation techniques are unlikely to be appropriate. 
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1 Introduction 



Systems that undergo first-order phase transitions as a function of temperature have 
been extensively discussed in the particle physics literature during the last decade. The 
initial motivation was the role of GUT-scale first-order phase transitions in driving 
an inflationary period of expansion in the very early universe [0. Lately the primary 
motivation has been to show that the observed baryon asymmetry of the universe can be 
produced through the baryon and lepton number anomalies, in concert with the out-of- 
equilibrium evolution of a first-order electroweak phase transition 0]. 

Most of the systems studied in these two contexts have been gauge theories where the 
first-order transition is radiatively induced [fluctuation-driven in the condensed-matter 
terminology we will employ in this paper). This means that the tree-level scalar potential 
has a form that implies a second-order phase transition, but when fluctuations in the fields 
are taken into account (at one-loop or higher) the potential indicates a first-order phase 
transition. If the transition is strongly first-order then the effective potential may be 
calculated by the standard loop expansion. However if it is weakly first-order (generically, 
when [Mh/M^Y ^ 1, where and are the Higgs and gauge boson masses) then 
the naive loop expansion breaks down, we must use the renormalization group (RG) 
to obtain reliable results 00. Ginsparg has investigated the Gell-Mann-Low RG 
flows of Yang-Mills-Higgs systems in the heavy Higgs regime, and found that the phase 
transition remains first order (specifically, there were no stable fixed points within the 
(4 — £:)-expansion) for reasonable scalar sectors. 

In this paper we propose to study the properties of weakly first-order phase transitions 
by the systematic use of the Wilson RG, together with an expansion in the number of 
spatial dimensions around four (the (4 — e)-expansion) 0. Following earlier work in 
the condensed matter literature, we will take the view that the appropriate quantity for 
the analysis of a first-order phase transition is the scale dependent ("coarse-grained") 
effective Hamiltonian defined by successively integrating out momentum modes. In the 
limit of infinite coarse-graining scale, this coincides with the effective action, and correctly 
describes the static properties of the system. However for a correct description of the 
dynamics we must coarse-grain only up to a physical scale, typically the correlation 
length. This approach has been advocated by Langer and coworkers 0, and by 
Kawasaki, Imaeda, and Gunton [||, and is also, of course, the strategy proposed by 
Wilson in his work on critical phenomena. We will explicitly calculate the coarse-grained 
effective potential for a toy model which possesses a fluctuation-induced flrst-order phase 
transition similar to that of the standard model. 

An important motivation for our RG study is that the current experimental lower 



1 



bound on the Higgs mass is around 60 GeV, so that the naive loop expansion for 
the electroweak phase transition is breaking down. Over the last few years there 
have been many attempts in the particle physics literature to remedy this situation 
[0]. They have mainly consisted of attempts to sum certain infinite subclasses of 
diagrams (for example "daisy" or "super-daisy" graphs, and "ring- improvement" ) and 
thus improve perturbation theory. We believe that the RG analysis is an advance over 
these diagrammatic methods for three reasons. 

Firstly, the RG analysis that we will discuss below generically leads to results for the 
properties of the phase transition (such as latent heats, bubble wall tensions etc.) that are 
highly non-analytic as a function of initial couplings (and e). This non-analytic behaviour 
is quite easily handled by the RG, but is very difficult to reproduce by diagrammatic 
methods. 

Second, as we will see in detail below, the typical flows of the effective couplings are 
such as to take them first towards, and then away from, an unstable fixed point. These 
unstable fixed points occur at values of the couplings of order 27r'^e. As e — > 1 (the 
physical value) this is far into the region of strong coupling where even the (non-Borel- 
resummed) renormalization group improved loop expansion is invalid. Diagrammatic 
methods obscure this problem rather than solving it. In this paper we will take e <^ 1 for 
most of our analysis, only setting e = 1 at the final stage. The motivation, and dangers 
[|], of this approach we will discuss at the end of the Introduction, and in Section 5. 

Thirdly, the Wilson RG flow of a theory down to a scale A is by definition the 
procedure which gives us the correct Hamiltonian for the dynamics of the modes of 
momentum \k\ < A, when fluctuations on scales \k\ > A are taken into account 
(integrated out). Thus, at a sensibly chosen scale, it is the relevant object for studying 
the dynamics of the phase transition. For example, the properties of the critical bubbles 
that nucleate the phase transition, such as the bubble wall tension, are determined by 
the effective Hamiltonian with coarse-graining length of order the correlation length 0. 
We show that the effective Hamiltonian at this coarse-graining scale calculahly differs 
from the usual effective action (which corresponds to an inflnite coarse-graining length). 

Rather than confront the full complication of the phase transition of a GUT or the 
electroweak model, in this paper we will study the following purely scalar system: 

1 11^ h h ^ 

^ = i^m) ■ id^v) + ^ E(^.)'- (1-1) 

Here Tp is an A^-component real scalar fleld, and Tp-Tp = Z^ili ^i^i- Inclusion of the second 
quartic coupling h2 breaks the original continuous 0{N) symmetry of the model down to 
the discrete hypercubic group in A^-dimensions. This system exhibits a close analog to the 
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weakly-first-order fluctuation driven phase transition that occurs in the standard model 
(and most other Yang-Mills-Higgs systems) for suitable Higgs mass. This will enable us to 
discuss some very important qualitative features of the physics of such phase transitions 
that have been missed by previous authors, without the inessential complication of gauge 
invariance. In the final section of this paper we will explain why our toy system is indeed 
a good theoretical laboratory for the general, gauged, case. 

There is one caveat that we should mention at this point. As one of has argued in an 
earlier publication the results of a 4 — e expansion can be unreliable; in particular, new 
fixed points of the RG can appear at e = 1 that are not accessible by such an expansion, 
rendering the transition second-order after all. However very useful information can still 
be extracted if the phase transition really is first-order, since its properties (such as 
the latent heat, critical bubble parameters etc.) can be calculated as an (asymptotic) 
expansion in e. Of course, the leading order calculation that we undertake in this 
work is not sufficient to accurately determine the properties at e = 1. However we 
are hopeful that higher-order calculations, in concert with asymptotic estimates of large- 
order e-perturbation theory, can supply accurate parameters in D = 3, as is the case for 
second-order phase transitions 0. Whether or not this succeeds, we have still extracted 
information about the limit e — > which other proposed calculations must be able to 
reproduce if they are correct. 

Before starting, we note that the RG has been previously applied to the cubic 
anisotropic model ( |l.lj ) by Rudnick [^], however he specialized to = 2, used 
classical statistical mechanics rather than finite temperature field theory, and used 
specific methods whose generalization to more complicated ( e.g. gauge) theories is not 
immediate; and also by Amit |TI| , who used the Gell-Mann-Low RG in which there is no 
finite coarse-graining scale. We have derived our evolution equations from Polchinski's 
general formulation using the perturbative expansion of Bonini, D'Attanasio, and 
Marchesini [|l^. This makes it relatively straightforward to push our results to higher 
order in the couplings, and even to include the evolution of derivative terms. 

Finally, we should mention that the philosophy of coarse-graining has also been 
advocated recently by Tetradis and Wetterich |T^, whose "average potential" is a 
smoothly cut-off version of the Wilson coarse-grained free energy (however, for a critique 



of the smooth cutoff see Ref. |^). They obtain accurate critical exponents for the 



second-order phase transition in an 0{N) model, but do not study first-order phase 
transitions. 
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2 The coarse-grained potential 



The Wilson RG is based on the idea of gradually integrating out degrees of freedom in 
a field theory by lowering the ultraviolet cutoff A, and at the same time varying the action 
functional in such a way that all Green functions of the remaining modes are invariant. 
We first describe how this is accomplished for zero temperature Euclidean quantum field 
theories, or equivalently classical statistical mechanical systems, and then generalize to 
finite temperature quantum field theories. 



2.1 The coarse-grained potential at zero temperature 

An elegant account of Wilson coarse-graining was given by Polchinski |2[ (see also Morris 



211]). Following his analysis we consider a simple scalar field theory defined by the 



functional integral 

Z[J] = y'[D0]exp|^[-i0(A;)(A;2 + m2)ir-i(A;VA')0(-A;) + J(A;)0(-A;)] +LaM 

= y"[(i0]exp5'A[(/)], 

(2.1) 

where K{k'^/A'^) is a smooth UV cutoff function that falls quickly to zero for k'^ > A^. 
The generating functional, Z, is independent of A (up to physically irrelevant field- 
independent changes of overall normalization) if the interaction Lagrangian L satisfies 

dm ^ _i / (27r)^ dK{kyA') f 6L 6L 6'L ] 

dA ^Jk{k^ + m^) dA \S(f){k)6(j){-k) 6(j){k)6(j){-k)} ^'^ 

together with J{k) = for /e^ > A^. The Green's functions of all the modes 0(A;) with 
^2 ^ ^2 then independent of A. No subtleties arise in the verification of these 
formulae since, with a physical cutoff, no divergences are present. The resultant coarse- 
grained effective action S'a[0] is the action that describes the dynamics of the degrees of 
freedom with momentum k < A. In the infinitely coarse-grained limit A — > it becomes 
the standard effective action. 

Note that the evolution equation ( |2.2| ) is an exact functional equation for the coarse- 
grained effective action. We wish to use this equation (actually a generahzation to finite 
temperature) to study the RG flow of our model system Eq. (|1.1|). It is, however, in an 
inconvenient form for our purposes - in particular the derivation of the "beta-functions" 
for the flow of the couplings is rather involved in this formulation. 

We will therefore employ the Wilson RG in a formulation originally constructed by 
Wegner and Houghton [|l^ in a statistical mechanical context, and further elaborated 



by NicoU, Chang, and Stanley [|l^. This formulation was then used to give the first 
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logically complete derivation of a coarse-grained effective action by Kawasaki, Imaeda, 
and Gunton (KIG) 0. (There are other discussions of this topic in the condensed- 



matter literature, see JT^ and [|T^], but overall there has been very little work on the 
derivation of the coarse-grained free-energy from a microscopic Hamiltonian.) Our choice 
of toy system Eq. ( |1.1| ) allows one major technical simplification in the passage from the 
Wilson-Polchinski equation ( |2.2| ) to the KIG form. Since our model involves only scalar 
fields with quartic interaction terms in the initial action, wavefunction renormalization 
only arises at two loops and above in the perturbative (coupling constant) expansion. 
This in turn implies that wavefunction renormalization is an O(e^) effect. However we 
will see that a calculation to 0{e) is sufficient to uncover the properties of the phase 
transition, and therefore, we can neglect the renormalization of the kinetic term to the 
order in which we work. Furthermore, all higher derivative terms generated as a result of 
the RG procedure are similarly higher order in the e-expansion. (Of course at very large 
momentum - in other words for field configurations with large spatial derivatives - the 
higher derivative terms can dominate the kinetic term. However for the configurations 
of interest, such as true vacuum bubbles, the gradients are such that this estimate is 
correct.) We can therefore ignore the renormalization of all momentum dependent terms 
- only the potential terms are renormalized to 0{e). This is a significant technical 
simplification relative to, say, scalar QED, or any other gauged system, where we would 
have to include wavefunction renormalization. We note in passing that it is possible to 
generalize our analysis to this case - this work will be reported in a future publication. 

In the case where we can ignore the flow of the derivative terms, the Wilson-Polchinski 
equation (P^ can be recast in the form: 



dA 2 



ff'V^ 



(2.3) 



This is our master equation for the flow of the potential of an A^-component scalar 
field theory in {D + 1) Euchdean dimensions. Here K(^d+i] = Sd/{'^'^)^^^^\ and Sd is 
the surface area of a unit D-sphere. This equation was first explicitly written down in 
Ref. 1^, so we will refer to it as the KIG-equation. In Appendix A we show that it can 
be derived from Ref 's. [|r^ perturbative analysis of Polchinski's exact RG flow equation. 
Intuitively, however, Eq. ( |2.3D seems very reasonable: it states that the contribution to 
the effective potential from integrating out the degrees of freedom in an infinitesimal 
shell A — (5A < < A in momentum space is just the standard loop result [|l9l for the 
effective potential, with the momentum integral restricted to that shell. 

The KIG equation as formulated above describes the flow of statistical mechanical 
systems, which correspond to zero temperature Euclidean quantum field theories. Given 
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the potential at some high cutoff Aq, it shows us how to ffow down to some scale of 
interest Ai and obtain Vai(0). Following KIG, we proceed in two stages: 

1. Expand Va(0) as a truncated power series in 0, including all the relevant operators 
(here defined to include the marginal operators as well). Recall that relevant 
operators are those whose couplings do not have large negative mass dimension 
for small e. We obtain coupled differential equations for the relevant couplings 
(only), which in the case of 0^ theory near four space-time dimensions are the 
mass-squared and the 0^ coupling(s) hi, which we solve to obtain m^(A) and 

2. Substitute the solutions of these equations into the RHS of the KIG equation, and 
integrate from Aq to Ai to obtain the higher order corrections to Vai, i.e. 0^ and 
higher. In general one must check for self-consistency — the "back-reaction" effect 
of these truncated couplings on the differential equations for the relevant couplings 
must be higher order in some small parameter. 

Note the standard 1-loop effective potential (with ultraviolet cutoff Aq and infra-red 
cutoff Ai) may be obtained by the crude approximation of skipping the first step, and 
substituting m(Ao) and 5'(Ao) directly into the RHS in the second step. 

In practical applications the theory would be specified by renormalization conditions 
on m? and the hi at a low momentum scale, so that at stage 1 we would choose values of 
m^(Ao) and /?.i(Ao) that give the right values when flowed down using Eq. ( ^.3| ), then at 
stage 2 we would calculate the remaining terms in the low-momentum potential. This is 
done because experimental measurements are typically made at energy scales low relative 
to the UV cutoff. In this paper, however, we are dealing with a toy model, and hence 
will express our results in terms of couplings defined at a high energy scale A^, defined 
in the next subsection. 

2.2 Coarse-grained potential at finite temperature 

The generalization of the KIG flow equation to finite-temperature becomes obvious if we 
recall the rules for constructing the finite-temperature (equilibrium) generating functional 
for a field theory: We compactify the Euclidean time direction on a circle of circumference 
13 = 1/T. This results in a discrete set of fourier modes Tp-,^[k.), with n E Z - the 
Matsubara modes. The (D + l)-momentum becomes k = (27mT, k), and integrals over 
are replaced by a sum over n. For example, the finite-temperature generating functional 
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for a simple theory, is, in momentum space: 



Z[J,T] 



I exp|-T ^ ^(-Kk' + (2™r)2 - m2)0„(k)0_„(-k) + J 
0„(k)0„,(k')0„"(k" 



y 

4! V 



k,k',k" 



-n—n'—n" 



-k - k' - k"l 



(2.4) 

It will be convenient to simplify this expression so that we have a conventionally 
normalized kinetic term - this is achieved by rescaling the field 0ri(k) by a factor of 
•s/T. With this choice, the quartic term is proportional to AT. This is convenient since a 
field theory in three dimensions has a quartic coupling of mass dimension one, which in 
this case is soaked up by T, so A stays dimensionless. 

In any case, the important point about Eq. ( |2.4] ) is that it is of the form of (infinitely) 
many coupled "species" of scalar field, indexed by n, living in one lower dimension. 
Therefore we can directly apply the KIG equation (|2.3|) to this system - actually the 
A^-component version - to find. 



dA 



{i,j),(n,m) 



In 



A jOn —m ~l~ r% r% 

d<^i,nd<^ 



j,m 



(2.5) 



as the appropriate generalization of the KIG equation to finite-temperature. Remember 
that in this equation the potential now includes a temperature-dependent effective 
(mass)^ term, (2'7rnT)^, for the n ^ modes. (We employ the tr(„ „^) notation as a 
shorthand for a sum over Matsubara modes where the discretized energy is conserved at 
each vertex.) 

There are two interesting limits of this formula. For A ^ 27rT, the Matsubara 
frequencies are so close that the "trace" over n approximates the original {D + 1)- 
dimensional integral, so Eq. (|2.3| ) is valid. For A ^ 27rT, the sum is dominated by 
the n = mode - the very massive n modes decouple leaving behind temperature- 
dependent corrections to the potential, and the system looks D-dimensional. In this limit 
Eq. (p^) thus simplifies to 



dA 



A^^'^hiij In 



A^6ij + 



d'VT 



(2.6) 



with the addition of a particular, temperature dependent, V^(^o) ^ boundary 
condition at the starting scale A^ < T. 

We must remember, however, that for a finite temperature field theory the 
renormalization conditions are still imposed on the zero temperature theory. The correct 
procedure is therefore to select some high scale Aq ^ T, and choose m^(Ao) and q^Aq) 
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such that the renormahzation conditions are obeyed when we flow down to low momentum 
in the zero temperature theory, i. e. using (p^). To calculate the finite temperature coarse- 
grained effective potential we should in principle run down from Aq using Eq. ( p.5|) , 
however the sum over Matsubara frequencies is difficult to perform analytically. 

Therefore in this paper we will instead arbitrarily choose some scale Aj somewhat 
below T, so that we can make the reasonable approximation of completely integrating 
out the massive n ^ modes. The running is then D-dimensional (in the physical 
case three-dimensional) below A^ and described by Eq. ( |2.(j| ). We will use the results of 
Ginsparg's perturbative calculation |^ of the effect of eliminating the non-zero modes to 
approximately calculate the resulting "boundary condition" potential at A^. Note that 
no infra-red difficulties arise in integrating out the n ^ modes as their propagators 
are cut off by large effective masses = (27rnT)^. Also note that we should not take 
At too much smaller than T or else the couplings of the effective Lagrangian will flow 
signiflcantly due to fluctuations in ^j Q(k), for k > At. Note that if greater accuracy is 
required, it is (numerically) straightforward to perform the correct procedure outlined 
above. 



3 Application to the cubic anisotropy model 

After this long preamble we are flnally ready to calculate the appropriate initial three- 
dimensional action at A^ for our toy system, and then give a general account of its 
fluctuation-driven flrst-order phase transition. Detailed analysis will be given in the next 
section. 

We are working at flnite temperature T, which enters the problem in two ways. Firstly, 
as we mentioned above, we rescale the flelds by a factor of VT so that all couplings in the 
action have their usual engineering dimension for a three-dimensional fleld theory, the 
four- dimensional couplings being modifled by factors of T. The mass term is unaffected, 
but the 0^ couplings are redeflned 

g,iA) = h,{A)T, (3.1) 

where hi are the dimensionless (4D) couplings. 

Second, the elimination of the n ^ modes shifts the values of the couplings /i^ and 



hi in Eq. (1J-), and also introduces higher dimension terms. However, only the correction 
to the value of /i^ is not suppressed by powers of the small couplings hi 0. For instance 
the shift in the 0"^ terms goes like hf - suppressed by one power of hi. For our model the 
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correction to the mass term is approximately 



2 , (iV + 2)/li+3/l2^2 



m\K,) = -/.^ + (3.2) 

where as we will see the first-order transition temperature is very close to the 
temperature where m?{T) = 0. 

Our project is therefore to analyse the phase transitions of the component scalar 
field theory defined by its truncated potential V^^ at (Eqs. (|3.3|) , ( p. 2D ) whose fiow 
into the IR is given by Eq. ( p.6|) . We will now follow stages 1 and 2 of the procedure 
suggested at the end of Sec. 2.1, first examining the fiow of V^^, then evaluating the full 
scale-dependent potential. 

The "truncated potential" , whose fiows we study, is 

vrm = \m\m ■ ^ + ■ + ^ Y:m\ m 

where we have replaced the coarse- graining scale A with a dimensionless variable /, 

A = exp(-/)A(0), (3.4) 

so / goes to infinity in the IR, and we choose our units of energy such that A(0) = Aj = 1. 
From the high temperature KIG equation ( p.6|) we obtain the fiow equations 

dm^ Kdb-^^ 1 



-({N + 2)g, + 3g2 



dl e~2' -|- m? 6 

dg2 Koe-'"' 



(2g^g2 + (3/2)(?2) • 



(3.6) 



dl (e-2' + m2)2 

A more useful set of variables is obtained by rescaling by the fioating cutoff, 

r{l) = m2(/)/A(/)2, 

A,(/) = g,{l)/A{lY^-^\ 

These are the variables in which a second order phase transition is described by a fixed 
point. Writing (4 — D) = e, the RG equations for {r, Ai, A2} are: 

I - 2K0 + ^i((iV+2)A, + 3A.) 
If = ^>'M-j^^,(>'i>'2 + {N + S)\V6) (3.7) 
W = -A.(0-^(2Aa.+ (3/2)A^). 
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Note that this is where we introduce the (4 — e)-expansion. The physical value is e = 1, 
but we will soon find it necessary to take e small in order to obtain any sensible results. 
For r <^ 1 we find four fixed points of (|3.7| ) (taking > 4 for the sake of simplicity - 
the < 4 case is entirely similar); they are (see Fig. 1), 

1. The doubly unstable Gaussian fixed point at Ai = 0, A2 = 0. 

2. The doubly stable Wilson-Fisher fixed point at Ai = 2e/{NKD), X2 = 2{N - 
A)e/i3NKD). 

3. The Ising fixed point at Ai = 0, A2 = 2e/{3K£)). 

4. The Heisenberg fixed point at Ai = Qe/{{N + 8)Kd), X2 = 0. 

The last two are both unstable in one direction in the {Ai, A2} subspace. 

As we will see below, first-order transitions occur when the renormalization group 
flows take the couplings out of the region where the quartic part of the potential is 
positive definite. There are two types of instability: 

1. Hypercubic diagonal instability. If Ai < then the potential is minimized for Tp 
in a diagonal direction, namely Tp = (±1, ±1, ±l)0/^/iV (of magnitude 0). The 
quartic part of the potential is then negative for Ai + A2/A^ < 0. 

2. Hypercubic axis instability. If A2 < then the potential is minimized for Ip along an 
axis, namely Tp = (0, 1, 0)0. The quartic part of the potential is now negative 
for Ai + A2 < 0. 

The unstable regions are indicated by shading in Fig. 1. 

Solving the fiow equations ( p.7|) with r ^ 1 yields the fiow pattern shown in Fig. 1. 
Note that there are regions of the coupling constant space that define sensible models ( 
i.e. that are bounded from below), that have fiows that inevitably take the couplings across 
the stability lines, and into regions without fixed points. For A^ > 4 these regions are the 
two wedges inside the stability region, but outside the quadrant Ai > 0, A2 > which is 
the domain of attraction of the Wilson-Fisher fixed point. Systems with {Ai(0), A2(0)} 
in these wedges will undergo fiuctuation-induced first-order phase transitions. Since the 
analysis of the phase transition for the two wedges is nearly identical, we will choose to 
concentrate on the lower wedge where A2 < 0, for which the RG fiows take Aj(/) across 
the stability line A2(0 = — Ai(/). In this case the phase transition is from a disordered 
phase to one with ordering along a hypercubic axis in ^-space. 

Also note that there exist regions with arbitrarily small initial couplings (Aj <^ e 
at At), that are first drawn towards the unstable fixed points, where the couplings are 
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of order e, before veering away and then crossing one of the stabihty hnes. Weakly 
first-order transitions correspond precisely to such values of the initial couplings where 
the amount of flow is large. Of course, there are also initial couplings that flow very 
little before crossing these lines, never becoming very large. These values correspond to 
strongly first-order transitions - for which it is not necessary to employ the RG to get 
a reasonable description of the transition, at least on infinite length scales. Ideally we 
would like to investigate the solution to the flow equations in the case e = 1. However 
the closed set of flow equations ( ^.7| ) was derived by truncation of an otherwise infinite 
set of coupled equations for the couplings in an arbitrary potential. Such a truncation 
is under control if the theory is weakly coupled for all scales A < At. For the weakly 
first-order case at D = 3 the couplings would quickly flow into a region of strong coupling 
where the loop-expansion and the closely associated truncation fails. This is the reason 
why a 4 — e expansion is necessary if we are to describe in a controlled way the weakly 
first-order case. 

In any event, for any given {r(0), Ai(0), A2(0)}, all three of order e or less, we can 
follow the RG flows using Eq. ( |3.7| ). Fixing Ai(0) and A2(0), we can then fine-tune r(0) 
so that the completely coarse-grained potential Ko(0) has two degenerate minima. This 
means we have chosen the parameters such that T is the temperature of a first-order 
phase transition. It turns out that we must tune r(0) to be very close to a critical surface 
^ci^i, A2), where ~ 0{e) (see Eq. (|4.1|)). As we evolve up in /, r(/) stays close to the 
critical surface, and hence the evolution follows the r <^ 1 trajectories of Fig. 1. At a 
value of / that we call we cross the stability line: 

Ai(/*) + A2(/*) = (3.8) 

We will see in Section 4 that in order to get the degenerate minima in V^o, we must have 
chosen r(0) such that r begins to deviate from Tc at h. (r(Z*) — rc{Xi{h), A2(/*)) ~ 0{e)). 
This deviation grows very fast, and at I = If (see Eq. (|4.27| )) r becomes of order 1: 

r{lf) = l + Oie), = + 1 ln(l/£) + 0(1) (3.9) 

At this point the anomalous scaling of r (the second term in the first equation of ( |3.7| )) 
is suppressed relative to the canonical scaling (the first term) by a factor of e. Since 
the anomalous term in the r flow equation has an factor of 1 -|- r(Z) in the denominator, 
it is now doomed: r grows like exp(2Z) for all / ~ //. The A's suffer the same fate - 
their anomalous scaling is suppressed by powers of 1 -|- r(/) too, so for / > //, {r, Ai, A2} 
flow canonically (or equivalently the physical parameters {m^ , gi, 92} are independent of 
/). This is illustrated in Fig. 1 by the straight line flow of the couplings {Ai, A2} in the 
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"unstable" regions. It does not matter that the A's will ultimately become large: the 
much faster growth of r suppresses one and higher-loop effects completely, and we can 
calculate reliably. These flows can then be substituted into the RHS of ( |2.6| ) to determine 
Vi for any /. This is the topic to which we now turn. 

If we take Tp to have magnitude along a hypercubic axis, the quadratic and quartic 
parts of the renormalized potential are 



V,'^{<P) = '^m'm' + j-,{g,{l)+g,{l)y 



(3.10) 



expressed in dimensionful variables. When we substitute this into the RHS of ( |2.6|) we 
get an integral expression for the higher- dimension terms that are generated as a result 
of the coarse-graining process: 



Vi{4>) = l^*'-(0) + AV^ {g^{s) + g^is))^^^ + {N - l)AVi 



(3.11) 



where, again we have specialized to an axis direction. Here AV is a functional, whose 
argument is a function of s over the range to I: 



AVi[w{s)] = ^Kd J^e-^^-'^'ds^^\n 



1 + 



Wis] 



Pis] 



^n£) , 1 

Pis) ^ 2 



W[S ) 

Ph) 



(3.12) 



Pis) = 



g 2s _^ m2(s) 



Note that we have subtracted off the quadratic and quartic parts of AVi, since these 
constitute the truncated potential, whose running we have already calculated. 

Iterating the KIG equation once again (before taking Tp to lie along the axis) we would 
find that the couplings of the induced higher- dimension terms change the RG equations 
for two and four-point terms. This change is perturbatively small when the four-point 
coupling is itself small. As we noted above, this is in general only the case when we work 
in a (4 — £:)-expansion with e small. 



4 The phase transition 

We now give the detailed proof of the qualitative results sketched above, and calculate 
the infinitely coarse-grained (/ = oo) effective potential near the phase transition, finding 
what value of m'^{0) yields degenerate minima for given gi{0), g2i0). 

4.1 The solution of the RG equations 

In view of the discussion at the end of Section 3.1, we only need solve the flow equations 
for I < If, for which the A's are always of order e. Ultimately we will show that the phase 
transition occurs at values of 771^(0) etc. such that (|3^) is satisfied. 
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We start with the flow of r. Defining 



rc(/) = -^((iV + 2)Ai(/) + 3A2(/)), (4.1) 

(note that Tc ~ 0{e)) we can rewrite the flow of r as 

1^2.-^. 0(e^ (4.2, 

As long as r ~ 0{6) we can ignore the (1 + r) denominator; When r ~ the whole 

second term is subleading. Thus 

nv 

— = 2r-2r^ + 0(er) yi<lf (4.3) 
dl 

Now dvc/dl ~ Ole"^), so to leading order in e, 

r{l) = rc{l)+ie^^ 
i.e. m2(/) = (4.4) 
t = e-2'/(l + 0(£)) 

Putting this into words, we have found that there is an unstable "critical surface" Vcil). 
If r starts off on the critical surface (t = 0) then it never leaves, its flow within the 
surface being determined by the flow of the A's. In order to see a subsidiary minimum in 
the potential we will have to choose i exceedingly small, so that If is large. Specifically, 
we must ensure that r does not grow to 0{1) until (just) after the couplings cross the 
stability line, i.e. at + 0{\ia{l/e). 

Our next task is to solve the flow of gi and g2, or equivalently Ai and A2. For I < l^, 
becomes 

dgi , /at , o\ 2 



= —^{Q9i92 + iN + 8)gi 



e'^f = -K,{2g,g, + {3/2)gl), (^.5) 



dx 

where x = exp{el) 



It is convenient, following Rudnick [|I0[, to define 



92 

[F is positive for the case we consider) so that 

dgi _ {N + 8)/3F^ - 2F 
dg^ ~ 3-4F ' 



F = -^, (4.6) 



(4.7) 
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and 

which is readily integrated to yield 

Now all we have to do is find F. By definition (f4.6|) and the RG equations (|3.5| ), 

— = ^(3F + {N~4)F') (4.10) 
ax 6e ^ ^ 

so by (pD, 

^ irD^2(0) i^(0)' (4.11) 

6e (3 + (iV - 4)F(0))" 
a = 3A^/(A^-4) 

Note that A is negative. Solving, we find the dependence of F on /: 

A{N - 4f {e'^ - 1) = [z(3 + {N -A)F{1)) - Z(3 + {N -4)F{0))} 

y2-a Qyl-a Qy-a (4.12) 



z{y) 



2 — a 1 — a a 



Armed with Eqs. (^4.9|) and ( [4.12|) we can now obtain gi{l) and (^2(0 ^ < ^* given 

^i(O) and (72(0). 



4.2 The effective potential near the phase transition 

We have managed to obtain the fiows in closed form. Stage 2 of our procedure now 
calls for us to substitute them in to Eq. (|3.11|) and send / ^ oo to get the effective 
potential, the form of which will tell us whether there is a first-order phase transition. 
In order to obtain understandable (and self-consistent) results we must again make the 
approximation of keeping only the leading order non-trivial terms in the e-expansion. 

We have previously argued (below Eq. (|3l9|) ) that r, Ai, A2 scale canonically for / > //, 
which means that m^, gi, and g2 can be taken to be constant. Thus the old non-RG- 
improved effective potential formula is valid in this region. It is therefore convenient to 
define the following variants of AV (see Eq. ( |3.12|) ) where the couplings are evaluated at 
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some particular / and do not flow: 

e-^^-'^'ds In 



P(s) 



1 + 



1 + 



w 



Pis) 
w 



P{s) 



w 



P{s) 



+1 



w 



Pis] 



(4.13) 



,-2s 



so that the potential becomes 



+ 



+ 92il))<P'] +iN- l)AVi^ [lgiil)<P' 
^Ui;'[\i9if + 92fW) + iN- l)A[///(|^?i;0^ 



(4.14) 



gif=giilf), etc. 



The first three terms represent the "ultraviolet" contribution, which requires integrating 
over the flowing couplings from / = to //. The last two terms are the "infrared" 
contribution, which are easy to evaluate, since nothing flows in the integrand. 

Actually, to the leading non-trivial order in e to which we work, one can also avoid 
the complications of integrating over flowing mass and couplings in the expression for 
AV (Eq. (|3.12| )), which is integrated from / = to //. To see this, first note that the 
flow of the mass in AV can be ignored, since it only enters through the propagator 
P{1) = exp(-2/) + m^{l) = exp(-2/)(l + 0{e)) + i, so we can set m'^{l) = m^ilf) = i 
to leading order in e, with i given in Eq. ( [4.4| ). Note that AV only supplies the (p^ and 
higher terms - the flow of m^, gi, g2 is explicitly taken into account in the lower order 
terms in V^"^. To see that the flows of gi and g2 can be ignored, look a typical term in 
AVy. 



D 



'2s 



ds 



in > 3) 



(4.15) 



where w{s) ~ g{s)(j)'^. Because t = exp(— 2//) this integral is dominated by the region 
close to //: the domain If — ^\ia{l/e) < I < If gives the dominant contribution to 
order e. In this relatively small range, the g^s hardly flow at all, since {1/ g)dg/dl ~ e, 
so we can take them to have their values at // throughout the whole range. Thus 
AV/^[t(;(s)] = AUV^ {w{lf)) to lowest order in e, and 



K»(0) 



Vi';i<p) 



+ AUry{\{g,f + ^2/)0^) + (iV - l)Af/,7(i^y0^ 
+ AUl^[\{g,f + g,fW) + {N- l)Af///(|<7i^0^ 



(4.16) 



We have now expressed the effective potential entirely in terms of simple integrals that 
depend on the relevant couplings at //. 
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Performing the integrals, we find that 



AUi{w) = AUl^{w) + AUl^^{w) 

= ^e^'-^^'hw - j^w' + fin + w)- fin) - iw^ + 2nw) ln(i + n)] 

fix) = x2(ln(x)-i), n=ril), 



w = e^'w, 



and defining a rescaled field 



(4.17) 
(4.18) 



we find that the potential is 



exp[(e - 4)//] + ^(Ay + A2/)0; + AUiA^Xif + A2/)0; 



+(iV-l)Af/,^(iAi;05) 



Ai/ = Ai(//), etc., 



(4.19) 



We are now almost finished. We can now see for what values of Aij,A2/ there is a 
phase transition at this temperature. We look for a subsidiary minimum at non-zero (pj, 
degenerate with the one at the origin. Anticipating the final result, we guess that it will 
occur for 

Ai/,A2/ ~ Oie), 
(A1/ + A2/) ~ OieHne) (4.20) 
0/ ~ Oil/e) 

Using this in ( |4.19| ) and keeping only terms of order (1/e^) Ine or 1/e^, we find 

Ko(0/) = + U^if + A2/)0/ + ^^^^^^^(|Av0j)^(ln(^Ai;0j) - 1) (4.21) 



which has the required degenerate subsidiary minimum for 
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[X^jKniN-1] 



(4.22) 



(which is of order 5^ Ine, as expected), at 



144 



XljKniN-1)' 



(4.23) 



(which is of order 1/e^ as expected). This demonstrates that indeed our model possesses 
a first-order phase transition of the type claimed. 
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We could at this stage perform the final task of re-expressing these values, and the 
mass r{lf) = 1, in terms of the initial renormalized parameters of the theory at ( 
i.e. / = 0) using Eqs. ( [4.9| ), ( |4.12| ), and ( ^4.4| ). Let us, however, delay this until our 
comments and conclusions in Section 5. 

It is also possible to express the infinitely coarse-grained free-energy in terms of the 
couplings evaluated at rather than If. The independence of the physical results (to 
leading order in e) enables us to check that we have not made an error in our calculation. 
To do this we need to fiow r^^., Ai/, A2/ up to First we must find which is defined 
by Ai(/*) + A2(/*) = 0. We therefore define 

u{l) = e-'^-'-'-^X, + X2) (4.24) 

and we find, by Eqs. ( p.7|) , (^^), to leading order in e, 



^ = -|A-.(7V-1)A;(0(^^ (4.25) 

We can integrate this equation from If to /*, since Ai varies by an amount subleading in 
e in this range. Under the correct assumption that exp(// — U) ~ 0{l/e), we obtain 

u{lf) = l,Kn{N - l)AL(i + ln[2e-2('/-'-)] + ■■■). (4.26) 

Hence 

which is much greater than 1 as anticipated. 

We can now go back to Eqs. ( glBD and ( gTzD and re-express the potential clS cl function 
of quantities defined at obtaining an expression like Eq. (|4.19|) with / — > *. Let us see 
for what mass and couplings we expect to see the degenerate subsidiary minimum. We 
flow the mass up to U by using (|4.4| ), ( [4.1|) and (|4.27|) 

rPT(/,) = ^ir^(iV- l)Ai,(ie-i/2 - 1). (4.28) 

The minimum will be at 



6e-V2 



Sh-i^){^yf = — . (4.29) 

Ai* 

We conclude, analogously to ( [4.201 ), that the interesting parameter range for the potential 



IS 



Xu = \2. ~ C(e), (4.30) 
0, ~ 0{e-'/^). 
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Keeping leading order in e we find the potential 



2r(/.)0, + - 



(4.31) 



This should be exactly the same as ( [4.21j ), and in fact it is. It has a degenerate subsidiary 
minimum at the values of r and predicted by ([4.28|) and ([4.29|) . Finally, note that 
our result is identical (to leading order in e and after suitable redefinitions) with that 
of Rudnick who previously analyzed a Landau-Ginzburg free energy of the form 
Eq. 



1.1) within a condensed- matter context. 



5 Comments and conclusions 

We now wish to make some comments concerning the formalism and calculation presented 
above. 

1) We have shown that it is possible to systematically calculate a scale- dependent 
coarse-grained free energy that describes a (fluctuation-induced) first-order phase 
transition. As we argued in the introduction, if one wishes to explicitly follow the 
dynamics of the phase transition (such as critical bubble nucleation and expansion) then 
one must coarse-grain not out to infinity (A 0), but to some appropriate length 
scale. We want to re-emphasize that this procedure should be familiar from the study of 
effective gauge theories - there one integrates out the heavy degrees of freedom to get the 
effective Lagrangian describing the dynamics of the light modes. In our case we might 
want, for instance, to accurately calculate bubble wall properties - we should therefore 
stop coarse-graining when we reach the bubble-wall thickness, and use the coarse-grained 
free energy at that length-scale. 

As an illustration of this procedure we can look at the (typical) case where this 
thickness is the same order as the correlation length. In terms of the variable ri used in 
Section 4 this is precisely when r/ = 1, i.e. / = If. We therefore need VJ^. However this 
is nothing more than the expression for the free-energy in Eq. ( f4.14|) without the extra 
infrared contributions 

VM) - ViM) = ^Uif{l{g,f + g^f)^') + {N- l)AU,f{lg,f<p'), (5.1) 
where we recall from ( [4.13 ) that 

roo , r ml 

(5.2) 



Af//^(w;) = \Kd I e'^^'-'^'dsln 



w 



In general, the difference between Voq{4>) and Vi^(0) (given for our toy model by Eq. (|5.1| )) 
can be quite significant. In our model it turns out that within the (4 — £)-expansion, and 
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for the orders of magnitude of 0, and the couphngs, of interest Eq. (|4.20|) , this difference 
is actually subleading by one power of e. Since we have not calculated our RG flows to 
higher order in e we therefore cannot give a reliable explicit expression for the values of 
the integrals in (|57T|). Of course, despite the fact that they are subleading in e, at £ = 1 
their contribution can be important. We expect that for the gauged case we describe 
below the difference between Voo(0) and VJ^. (0) will not be subleading in e. 

2) We now want to return to the one loose-end of Section 4 - namely the re-expression 
of the properties of the potential in terms of the initial variables at / = 0. Using the 
solutions, Eqs. ([4. 91), ( [4.12| ), and ([4.4]), we derived in Section 4.1 to the RG equations this 
is a straightforward task. However the expressions one gets in doing this are quite lengthy. 
A simplification occurs if we take the initial ratio of the couplings Fq = —gi{0)/g2{0) ^ 1. 
For example, in this limit the value of the field at the new (true vacuum) minimum is 
given by 



^ N{N + 5) + 3 £ 



where 



_ f NjN + 2)iN + 8) gMKn (N - -^^^^ 

-{ iv(iv + 5) + 3 e^liv^j / ^° • ^^-^^ 

(/* is the value of the flow parameter at which the system crosses the stability line 
Qii}*) + g2ij'*) = 0). As advertized in the Introduction these expressions are quite non- 
analytic as a function of the initial couplings and e. The expression for 0^ Eq. ( p.3| ) also 
exemplifies the connection between weakly first-order transitions and significant flow in 
/ that we mentioned in Section 3. This is because when is large the value of 0^ in the 
new minimum (at the transition temperature) given by ( ^.3|) is small compared to its 
zero-temperature value. 

3) In the introduction we motivated the study of our toy scalar model by stating 
that it shared many features in common with the phase transition occurring in gauged 
systems. One way of understanding why this is true is to look at the RG flow equations 
(once again in (4 — £:)-spatial dimensions) for the couplings of the prototypical gauged 
system - the iV-component Abelian Higgs model. Again we have two couplings that 
are marginal near four- dimensions, the quartic scalar self-coupling h{l), and the gauge 
coupling e^(Z). Their RG equations are 

rfe^ , NKd a 

^ f5 5) 

dh , . {N + A)Ko ^, , ^ ■ ' 
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where t] = —'iK£,e^/{l + r) is the anomalous dimension of the scalar field induced by 
wavefunction renormalization, and r{l) is the (dimensionless) scalar mass parameter with 
a RG equation we shall not display. (As in Section 3 we have implicitly included the 
appropriate factor of temperature T induced by the rescaling of the fields, as well as then 
dividing out by A(/)^ so as to get dimensionless couplings.) For any reasonable number 
of complex scalar fields N < 366 these equations have only two fixed points - the doubly 
unstable Gaussian fixed point at the origin, and a fixed point at 

= 0, (5.6) 



which is unstable in the direction of the gauge coupling. The qualitative pattern of the 
fiows in the {/i, e^} plane is illustrated in Fig. 2. 

The important point is that there are well-defined models (indeed the whole of the 

> 0, /i > quadrant!) whose fiows inevitably take the system outside the naive 
stability region (the same > 0, /i > quadrant). Just as in the scalar case we have 
considered above, this implies that the system undergoes a fluctuation induced (finite- 
temperature) phase transition. Furthermore, the work of Ginsparg shows that the 
standard model and almost all other gauge systems of interest fall into this class. In a 
future publication we will show how the properties of these first-order phase transitions 
can be analyzed by the application of the Wilson RG and the (4 — £:)-expansion. 

4) In a recent very interesting work, Shaposhnikov |20| argued that to fully understand 
the finite-temperature phase transition that occurs in gauge theories, it is necessary to 
include non-perturbative contributions to the free-energy. This, he argued, is because 
the effective three-dimensional theory of the zero Matsubara modes is a confining theory 
in the far infrared - at least in the non-Abelian case. It was also argued that the 
non-perturbative effects can make the transition in, say, the minimal standard model 
considerably stronger for a given Higgs mass than a perturbative analysis would indicate. 
It is interesting to see what our approach says about this general question. 

We have already noted that the effective (4 — £)-dimensional theories that we must 
consider if we want to study the phase transition have the property that as we flow into 
the infra-red the effective couplings flrst flow towards the unstable flxed point, which 
for e = 1 is well into the strong-coupling regime. The same phenomenon occurs in 
gauge theories, as exemplifled by the RG flows of the Abelian Higgs model illustrated in 
Fig. 2. One is therefore tempted to say that this is almost a proof of the correctness of 
Shaposhnikov's suggestion at £ = 1. All the theories are in their strong-coupling regime 
and non-perturbative effects are important. However, this does not necessarily mean 
that we cannot calculate. 
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To see this, consider the flows in Fig. 1 that take us towards the stable Wilson-Fisher 
fixed point at Ai = 2e/ (NKd), A2 = 2(iV — 4)e/ {3NKd), corresponding to a second-order 
transition. For e = 1 this is also far into the strong- coupling region, so surely we must 
consider non-perturbative effects if we are to get accurate predictions for the properties 
of this transition (such as the critical exponents). In fact, the experimentally observed 
critical exponents are amazingly well described by an expansion up to only O(e^). Even 
more accurate results can be obtained by Lipatov techniques for estimating large order 
behavior of the e-expansion together with Borel resummation. These techniques are 
closely related to the usual semi-classical methods employed in dealing with weakly non- 
perturbative effects. Thus we might say that for this second-order case it seems to be 
that the non-perturbative effects that are present are fairly benign - and can be handled 
within more sophisticated versions of the e-expansion. 

What about the first-order case - especially for gauge theories? Here the situation 
is less clear, mainly as a consequence of the relatively limited amount of work that been 
performed so far. There has been for example, to our knowledge, no extension of the 
Lipatov techniques to the first-order case. One obvious difference between the first and 
second-order cases is that in terms of the dimensionless variables, Xi{l), the flows take 
one out to "infinity" in coupling constant space, rather than to a fixed point. This is 
actually benign in the scalar case since soon after we cross the stability line we reach 
a point, //, where the increase in the scalar mass cuts off any further non-canonical 
scaling of the couplings, so the Aj(/) move along a straight line trajectory projecting 
out from the origin (see Fig. 1) - in other words the dimensionful couplings are to a 
good approximation constant after If. However, in the gauge case the flow of the gauge 
coupling doesn't seem to be cut off in quite such a simple way. This is consequence of 
the lack of any factors of 1/(1 + r(/))"' appearing in the RG equation for in ( ^.5|) . Non- 
Abelian theories also present new difficulties. We are at present considering these issues. 
In any case it is certainly true that the RG techniques we have utilized in this paper 
lead to considerably greater insight into this question than more standard diagrammatic 
resummation techniques. 
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Appendix A Derivation of the evolution equation 



In this appendix we will derive the KIG evolution equation ( |2.3| ) from Polchinski's 
evolution equation ( |2.2| ), for the case of a single scalar field in 4 dimensions. The 
generalization to multiple fields and arbitrary dimension is straightforward. We first 
show that, for momenta below the cutoff A, Polchinski's interaction Lagrangian 
is the same as the IPI generating functional T\[(f)\ defined by Bonini, D'Attanasio, and 
Marchesini (EDM) |jl3| in their perturbative analysis of Polchinski's equation. We then 
use their evolution equation for T/<^[(f)] to derive the KIG evolution equation for the 
effective potential. 

To clarify matters we will take the momentum cutoff at A to be sharp, and use 
superscripts > and < to distinguish degrees of freedom above and below the cutoff 



respectively. (For a rigorous analysis see Morris's very clear exposition |21], in which 
the sharp cutoff limit is taken at the end.) The propagator has two domains. 



Infrared (IR): D^{k) 
Ultraviolet (UV): Dy{k) 



K{kyA^)/{k^ + m^), 
[1~ K{e/A^)]/{e + m^) 



(A.l) 



These add up to the full propagator D{k). For most purposes it is natural to take K{x) 
to be a step function that is 1 for x < 1 and zero otherwise. To make sense of EDM's 
definitions, however, we will have to allow the UV propagator to have some infinitesimal 
value in the IR, and take the limit as it goes to zero. 

The full generating functional of the theory can be split into IR and UV portions: 



Z[J] = Z[J<,.P 



[D0]exp / \-U{k)D-\k)<P{~k) + Jik)^{-k)] + L^,,, 



[D(p<] exp 



U<{k)D-\k)<p<{-k) + J<{k)<p<{-k) + La[0<, J>] 



(A.2) 

Comparing with Eq. ( |2.1|) we see that Polchinski's interaction Lagrangian La is the 
functional integral over UV fields, with a given IR background field 0^ and in the presence 
of a UV source J^, 



exp L/ 



[D0>] exp 



\(P>ik)D-\k)(P>i-k) + J>{k)(j)>i-k) 



+ L 



bare I 



(A.3) 
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Eq. ( |A.2| ) tells us that we should think of La[0^, 0] as the effective action "felt" by the IR 
fields. We now want to show that it is the same as EDM's effective action (for momenta 
below A) 

T4^P<,J> = 0]= [ j<ik)^<ik)-W[j<,Ol with V^< = ^, (A.4) 

Jk<A OJ^ 

which is the Legendre transform of their connected generating functional with UV 
propagator, 

expWAir, J>] = [ [D0]exp( / \-U{k)D-\km-k)+j<{k)<f)<{-k) + J> {k)<j)> {-k) 



+ Lbare[0<,0>] 

(A.5) 

A peculiar (but essential) aspect of this definition is that the functional integral is over IR 
as well as UV fields, even though the propagator is cut off in the IR. Moreover, there is an 
IR source coupled to the IR field. This only makes sense if we give D> some infinitesimal 
value e/lk"^ + m^) in the IR, and perform the IR part of the functional integral before 
sending e ^ 0. Dy^{k < A) is then very large, and in combination with the j^(f)^ term 
gives a 5-function, requiring (p^lk) = tp^{k) = D~^{k)j^{k). The result is 

exp-rA[V^<,0] =expf-i / ^<{k)D-\k)^<{-k)) 

X j[D(t)>]exp)^jy\<P>{k)D-\k)<P>{-k)+\ -Lbarc[^<,0>]|. 

We see that rA[V'^] is indeed the same as Polchinski's effective action La[</>^] Eq. ( [A. 3D , 
except that it has an additional inverse propagator term in front. As we send e — 0, 
eliminating EDM's infinitesimal IR component of D>, this term goes to infinity {D^{k < 
A) ~ e). This is not a disaster — it just refiects the fact that we in effect gave the IR 
modes infinite mass in order to treat them as a background. It only affects the 2-point 
function, which EDM write as 

(A;) = D-^(A;) + SA(fc), (A.7) 
and they obtain an evolution equation for the e- independent part Sa. We therefore write 

Va(0) = (m2 + EA)02 + FA(0), 

OO 1 



n=2(2n)! ^ 



where the are zero-momentum n-point IPI functions. 
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We now begin the second stage of our calculation, showing that EDM's evolution 
equation for Va reproduces the KIG equation Eq. (|2.3|) , in the approximation where 
momentum-dependence in the action is ignored. Expanding Va in the same way as Fa, 
so that ^ = m2 + Ea but forn > 1, = ri^"\ EDM find 

where we have written D\{q) = D~^{q), and 



ff +^)(g) = rf - V ]^^^^ j^, (A.IO) 

^ ^ /^A'(?) + SA(g) ^ ^ 

so that 

'^^^ Af=l I /i=l l2=l In-i=1 



(2/i)!---(2Z^r)! 



(A.ll) 



^J- A ■ ■ ■ J- A 



The momentum integral is contained in In, 

Assuming that the step function K{q/A) (see ( [A.l|) ) varies very quickly near g = A, we 
can evaluate it, 

In{A) = ^{A' + m' + Sa)-^. (A.13) 

On the RHS of Eq. ( [A.ll ), the quantity in braces contains — 1 summations. If 
= 1 then it evaluates to the quantity in square brackets, with 1^ = n. Otherwise 
the summations are evaluated in the usual way, with 1^ = n — li — ■ ■ ■ — l^-i- Summing 
over 77,, 

QT/ oo _( _'\\N ( CO p2;i+2 oo p2Zjv+2 ~v 

(A.14) 

so by (^), 

A^^iK.AMn(l + -im^). (A.15) 
OA 2 \^ A2 + m2 + SA/ 

Ignoring terms that are independent of (f), 

A^^ = iir4A^ln(A^ + V:('(0)), (A.16) 

which was to be proved. 
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Figure Captions 

Figure 1. The general form of the RG flows in the {Ai,A2} subspace are displayed, 
along with the four fixed points, for the case > 4. Shading indicates the "unstable" 
regions where the quartic part of the potential is no longer positive definite. The values 
I = h and I = If of the flow parameter are where, respectively, a stability line is crossed, 
and r = 0{1) so that the subsequent flow is canonical ({Ai(/), A2(/)} flow along a straight 
line projecting through the origin). 

Figure 2. The general form of the RG flows in the {/i, e^} subspace for the A^- 
component Abelian Higgs model for A^ < 366, and in (4 — £:)-dimensions. Both fixed 
points are unstable. The physical "unstable" region is /i < 0. 
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